Absence of orbital-selective Mott transition in Ca 2 _ a; Sr :! ;Ru04 
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Quasi-particle spectra of the layer perovskite Sr2Ru04 are calculated within Dynamical Mean 
Field Theory for increasing values of the on-site Coulomb energy U. At small U the planar geometry 
splits the tig bands near Ef into a wide, two-dimensional d xy band and two narrow, nearly one- 



dimensional d x 



bands. At larger U, however, the spectral distribution of these states exhibit 



similar correlation features, suggesting a common metal-insulator transition for all t^g bands at the 
same critical U. 



Sr 2 Ru04 has attracted considerable interest during the 
recent years because of a variety of fascinating prop- 
erties. It exhibits unconventional p-wave superconduc- 
tivity P, and is in fact the only known layered per- 
ovskite without copper that becomes superconducting in 
the absence of doping. The ground state is non-magnetic 
whereas the volume perovskite SrRu03 is a strong fer- 
romagnet. Substitution of Sr with Ca makes Sr 2 Ru04 
undergo a transition to an antifcrromagnetic Mott insu- 
lator In contrast, doping with La leads to a fer- 
romagnetic instability accompanied by a breakdown of 
Fermi-liquid behavior [{|. Spin fluctuations in Sr2Ru04 
are believed to have both ferromagnetic and antiferro- 
magnetic components 0. Moreover, inelastic neutron 
scattering data indicate a soft phonon mode associated 
with in-plane rotations of oxygen octahedra |jj • Clearly 
Sr2Ru04 is close to structural and magnetic instabilities. 
Substitution of Sr by the smaller Ca ions also causes a 
crystallographic distortion consisting of tilting and rota- 
tion of oxygen octahedra. These modifications reduce the 
effective hopping between Ru Ad orbitals via O 2p levels. 
The resulting narrowing of the Ad bands is believed to be 
the driving mechanism for the metal-insulator transition 
observed in Ca 2 -xSr x Ru0 4 for x = 0.2 0, HQ. 

The layer perovskite Sr2Ru04 can be regarded as a 
prototype of a strongly correlated, highly anisotropic 
multi-orbital transition metal oxide with coexisting wide 
and narrow bands near the Fermi level. Due to the pla- 
nar crystal structure the partially filled £ 2s bands sepa- 
rate into a two-dimensional d xy band with a van Hove 



to within approximately 10 meV of Ep. Thus, doping 
with a few % La moves the Fermi level through the d xy 
singularity, thereby causing deviations from Fermi liquid 
behavior y|. 

At the surface of Sr 2 Ru04, a y/2 x \pl lattice recon- 
struction driven by a soft phonon mode involving rota- 
tions of O octahedra reduces the effective d — d hopping 
|15| . As a result of the slightly narrower width of the d xy 
band in the first layer the van Hove singularity is shifted 
below Ep, converting the 7 sheet at the Fermi surface 
from electron-like to hole- like. Thus, photocmission spec- 
tra reveal a superposition of bulk bands with E v h s > Ep 
and surface bands with E v h s < Ep 0] • The bulk bands 
are then consistent with the Fermi surface derived from 
de Haas- van Alphen measurements |Tt| . 

To investigate the metal-insulator transition induced 
by iso-electronic substitution of Sr with Ca, Anisimov 
et al. recently performed DMFT calculations for 
Sr 2 Ru04 for increasing values of the local Coulomb re- 
pulsion U. While in reality the transition is driven by 
a distortion of O octahedra, giving rise to a narrowing 
of the Ru Ad band without any significant increase in 
U, these authors argued that a qualitative understand- 
ing of the Mott transition can be achieved by increasing 
U and ignoring the modifications of the density of states 
due to the structural transition. This picture is very in- 
teresting from a conceptual point of view since in the 
anisotropic £ 2g configuration U is comparable at first to 



W x 



z,yz 



1.3 eV and subsequently to W xy « 3.5 eV. A re- 



markable result of these calculations were two successive 
singularity (vHs) just above Ep |9j, and two nearly one- Mott transitions for the d XZtVZ and d xy bands, implying 



dimensional d xz>yz bands. In the single-particle picture 
these bands do not hybridize for symmetry reasons. Lo- 
cal electron-electron interactions in the Ru Ad shell, how- 
ever, are strong and according to angle-resolved photoc- 
mission spectra the on-site Coulomb energy lies be- 
tween the single-particle band widths of the d xz . yz and 
d xy bands: W xz ,yz <U < W xy 0]. Thus, although cor- 
relations primarily affect the narrow d xz . yz bands, multi- 
orbital interactions also influence the wider d xy bands, 
in particular its singularity above Ep. Recent quasi- 
particle calculations [ll| based on the Quantum Monte 
CarlofQMC) Dynamical Mean Field Theory (DMFT) 
(T2I lill YIA\ showed that on-site correlations give rise to 
a charge transfer from the narrow bands to the wider 
d xy band and to a red-shift of the van Hove singularity 



■xz,yz 

an intermediate region of U or, equivalently, of Ca con- 
centration, where the narrow bands exhibit an excitation 
gap while the wide band is still metallic. Only at still 
larger U a gap appeared for the entire t 2g complex. 

To study this unusual electronic structure in greater 
detail we have extended our previous DMFT calculations 
for Sr 2 Ru0 4 E3 to larger values of the on-site Coulomb 
energy. Surprisingly, we find a behavior that differs fun- 
damentally from the one reported in Ref. 0. For small 
U, the d xy and d XZtVZ quasi-particle spectra indeed have 
quite different shapes owing to their predominantly two- 
versus one-dimensional electronic structure, respectively. 
As U increases, however, these spectral distributions be- 
gin to resemble one another, with only minor differences 
in the weight and shape of the coherent peak near Ep 
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FIG. 1: (a) Tight-binding fit to LDA t 2g bands of Sr 2 Ru0 4 
(Ep = 0); (b) Density of states Pi(u>) of d xy bands (solid 
curve) and of d XZiVZ bands (dashed curve). 



and of the lower and upper Hubbard bands. Our results 
therefore suggest that for sufficiently large U the corre- 
lated electronic structure within the ti g shell dominates 
the one-electron hopping. As a consequence, one com- 
mon metal-insulator transition exists at an intermediate 
critical value of U between those of the isolated d XZ:VZ 
and d xy bands. The reason for the discrepancy between 
the two approaches is not clear at present, but could be 
related to the fact that in Ref. Q the quantum impurity 
problem was treated within the non-crossing approxima- 
tion (NCA) [H m whereas we use the more accurate 
Quantum Monte Carlo method [l3j . 

Fig. 1 shows the tight-binding fit to the LDA band 
structure of Sr2Ru04 Q in the vicinity of Ep and the 
density of states. In the bulk the d xy vHs is about 50 meV 
above Ep. All three bands are approximately 2/3 filled. 
The non-symmetric shape of the two-dimensional p xy (uj) 
is caused by second-neighbor Ru-Ru hopping terms. The 
d xz ,yz densities are dominated by one-dimensional hop- 
ping along the Ru rows. 



To take account of local Coulomb interactions we cal- 
culate the quasi-particle spectra within the multi-orbital 
QMC-DMFT approach [ill HJ. Since the d xy and 
d xz ,y Z bands do not hybridize, the self-energy is diag- 
onal in orbital space. In the DMFT the elements Ej 
(i = xy, xz, yz) are functionals of the bath Green's func- 
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tions G- 



G i + Ej, where the local Gi is given by 



and n is the chemical potential. Most QMC calcula- 
tions were done for (3 = 8 with 64 time slices and sev- 
eral runs using 2xl0 5 sweeps. In the low temperature 
calculations for (3 — 40 128 time slices were used, with 
up to 5x10 sweeps. The quasi-particle density of states 
Ni{ui) = —lmGi(uj)/TT is obtained via maximum entropy 
reconstruction |2l| . 

As shown in Ref. ^l| > the angle-resolved photoemis- 
sion spectra of Sr2Ru04 can qualitatively be represented 
by quasi-particle distributions for local Coulomb and ex- 
change energies U — 1.2 eV, J = 0.2 eV. The d xy and 
d xz , y z spectra then look like moderately deformed ver- 
sions of their respective single-particle densities of states. 
For instance, the lower van Hove singularity of the d xz ^ yz 
bands near —1 eV is broadened and shifted to about 
—0.5 eV. Larger Coulomb energies would yield too small 
binding energies for this spectral feature. This result is 
consistent with the fact that on-site Coulomb energies 
for 4d transition metals are smaller than for 3d metals. 
The d xy vHs is also shifted and lies only about 10 meV 
above Ep. Both bands show very low spectral weight 
in the energy region of the Hubbard peaks. The effec- 
tive masses calculated within the QMC-DMFT for these 
Coulomb energies also agree with experiment. 

We now increase the on-site Coulomb energy in or- 
der to explore the metal-insulator transition for this 
anisotropic multiband system. To illustrate the effect 
of the different widths of the t2 g subbands on the quasi- 
particle spectra we show first in Fig. 2 the results for 
hypothetical three-fold degenerate bands consisting ei- 
ther of d xy or d xz ^ yz character. The total filling is 4 as 
in actual Sr 2 Ru04. For U = 3.0 eV the wide d xy band 
is dominated by the strong quasi-particle peak at Ep. 
Since U < W xy it is only moderately affected by correla- 
tions and shows only weak shoulders in the range of the 
Hubbard bands. On the other hand, since U >> W XZtVZ 
the narrow band derived from the d xz ^ yz density of states 
exhibits clear signs of a Mott transition. The upper and 
lower Hubbard peaks are the dominant spectral features. 
Because of the finite temperature used in the QMC- 
DMFT calculation, the remaining quasi-particle peak at 
Ep obscures the gap between the Hubbard bands. Nev- 
ertheless, according to the temperature dependence of 
similar features found in one-band systems 0, |53 the 
value of U used in the spectra shown in Fig. 2 should 
be very close to the critical value for a metal-insulator 
transition within the d xz . yz bands. 

The quasi-particle spectra for the actual ti g bands of 
Sr2Ru04 are shown in Fig. 3 (a) for the same on-site en- 
ergies. Evidently, as as result of the strong inter-orbital 
Coulomb and exchange interactions both d xy and d xz ^ yz 
spectral distributions have nearly lost their original den- 
sity of states character and look remarkably similar: with 
coherent peaks at Ep and upper and lower Hubbard 
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FIG. 2: Quasi-particle density of states Ni(u) of Sr^RuCU 
derived from DMFT for T = 1450 K, assuming indepen- 
dent three-fold degenerate (a) d xy and (b) d X z,yz bands (solid 
curves; see text). U = 3.0 eV, J = 0.2 eV; dashed curves: 
bare densities. 



FIG. 3: Quasi-particle density of states Ni(ui) of Sr2Ru04 
derived from DMFT for T = 1450 K. (a) U = 3.0 eV, J = 
0.2 eV; (b) U = 4.0 eV, J = 0.7 eV. Solid curves: d xy states; 
dashed curves: d xz ^ yz states; dotted and dash-dotted curves: 
corresponding bare densities. 



bands of about the same intensity, position and shape. 
These spectra suggest a common Mott transition for the 
t2 g bands at a value of U only slightly larger than 3 eV. 
We point out that the local correlations for U = 3 eV 
induce a charge transfer of about 0.05 electrons from the 
narrow d xz _ yz bands to the d xy band. This would im- 
ply that the d xy van Hove singularity has shifted below 
Ep and that the 7 sheet of the Fermi surface is hole- 
like rather than electron-like as observed in de Haas-van 
Alphen measurements EJ- This supports our previous 
choice of a smaller U for the actual Si'2Ru04 t2 bands 

El. 

Fig. 3 (b) illustrates the effect of increasing the ex- 
change energy to J = 0.7 eV. In order to keep the average 
Coulomb energy U unchanged U is increased to 4 eV. 
(For a tig complex U coincides with the inter-orbital 
Coulomb energy U' = U - 2 J @]. Thus, U = 2.6 eV in 
Figs. 3(a) and (b).) Although the d xy and d XZiVZ spec- 
tra are not quite as similar as in Fig. 3(a), they have 
coherent and incoherent peaks of about the same inten- 
sity. Only the positions of these features are shifted. Ex- 
trapolation to lower temperatures does not seem to be 
quite as straightforward as in the previous case. Nev- 
ertheless, qualitatively both d xy and d xz , yz spectra are 
equally correlated and therefore also suggest a common 
Mott transition. 

Judging from the results discussed so far for U in the 
critical region it does not seem likely that upon lower- 



ing the temperature the d xy and d xz ^ yz spectra would 
begin to lose their similarity again and diverge towards 
two separate Mott transitions. To test this behavior 
we have performed multiband QMC-DMFT calculations 
at about 300 K which are computationally much more 
costly. Fig. 4 shows the d xy and d XZjVZ quasi-particle dis- 
tributions for U = 3.0 eV, J = 0.2 eV. The peak near E F 
exhibits a fine structure which is the remnant of the up- 
per and lower edge singularities of the d xz%yz bands. The 
d xy spectrum also reveals this structure as a result of the 
strong orbital interactions. The important point is that 
both d xy and d xz ,yz spectra are similarly correlated: the 
coherent peak has about the same weight and the upper 
Hubbard bands are almost identical. The main differ- 
ence is that the lower Hubbard d xz ^ yz satellite is split off 
from the main peak by a gap while the d xy satellite lies 
in the lower tail of the d xy density of states and therefore 
forms a continuum with the main peak. Leaving aside 
the uncertainties caused by the maximum entropy recon- 
struction, such a gap does not imply that the spectral 
weight at the Fermi level vanishes at a different U CT it or 
T cr it than in the case of the d xy states. On the contrary, 
the comparable strength of the coherent peaks indicates 
that the d xy and d xz . yz bands undergo a common Mott 
transition. 

The above results differ qualitatively from the ones 
obtained by Anisimov et al. 8]. The NCA in principle 
should be applicable also to multiband materials. On the 
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FIG. 4: Quasi-particle density of states Ni(u) of S^RuCU 
derived from DMFT for T = 300 K, U = 3.0 eV, J = 0.2 eV. 
Solid curves: (a) d xy states; (b) d XZ:VZ states; dashed curves: 
bare densities. 



other hand, even in one-band cases it is known to be unre- 



liable in certain situations and even violate Fermi-liquid 
behavior 0, 0, |2£j . It would be desirable to perform 
numerical renormalization group calculations for this sys- 
tem at very low temperatures. Unfortunately, multiband 
calculations within this scheme are computationally not 
yet within reach. 

In summary, we have used the QMC-DMFT to cal- 
culate quasi-particle spectra of Sr2Ru04 in the range 
of on-site electron-electron interactions appropriate for 
a possible metal-insulator transition. Qualitatively, in- 
creasing U might be a way to simulate the reduced ti g 
band width caused by distortions of O octahedra when Sr 
is replaced by Ca. At low Coulomb energies, the spectral 
distributions of the d xy and d XZiyz bands differ strongly 
because of the planar geometry of this perovskite mate- 
rial. Despite this anisotropy, at sufficiently large U these 
spectra resemble one another in the sense that they are 
equally strongly affected by local correlations. In par- 
ticular, they exhibit similar upper and lower Hubbard 
bands and a coherent peak of similar strength. These 
results suggest that there is a common Mott transition 
for the tig complex, in contrast to the squential orbital- 
selective Mott transitions found in Ref. [8j . Thus, in the 
critical region of U local Coulomb interactions appear to 
dominate the anisotropic one-electron hopping. 
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